Quasicondensation in 2D Interacting Bose Gas: Quantum Monte Carlo Study 
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We present a detailed Monte Carlo study of correlations in 
an interacting two-dimensional Bose gas. The data for one- 
particle density matrix in coordinate representation are com- 
pared to the results for the local many-particle density cor- 
relators, which are responsible for the recombination rate in 
real experiments. We found that the appearance and growth 
of quasicondensate fluctuations changes local correlations be- 
tween the particles well before the Kosterlitz-Thouless transi- 
tion point. The amplitude of the l/m!-effect is very sensitive 
even to a rather moderate interaction, and is considerably 
smaller than its limiting value. 

PACS numbers: 03.75.Fi, 05.30.Jp, 67.65. +z 

The discovery of Bosc-Einstcin condensation (BEC) in 
ultra-cold dilute gases has opened a unique possibility 
for the study of quantum correlations in these systems. 
For one thing, it was observed that inelastic processes are 
suppressed in the presence of BEC jjj , in agreement with 
theoretical predictions Q . In the 3D gas at low temper- 
ature T <C T c (T c is the BEC transition temperature), 
the m-particle correlator 

J rr„ l = n- m ((vl/t(0,0)) m (*(0,0)r) , (1) 

is reduced by a factor of 1/m! as compared to its value 
at T > T c Here "J is the Bose field operator, and n is 
the gas density. This result is exact for the ideal gas, and 
interaction corrections are small in the gas parameter. In 
particular, the rate of three-body recombination (which 
is proportional to K3) must drop by ~ 6 times 0], as was 
indeed measured in Ref. 0. 

In systems with reduced dimensionality (e.g., in d = 2) 
at any finite T the condensate density is zero, and the 
question of l/m!-effect is not at all obvious. On another 
hand, below Kosterlitz-Thouless (K-T) temperature T c 
2D system becomes superfluid due to long-range phase 
coherence. In the superfluid state one finds that phase 
and density correlation lengths, R c and r c , may have 
different scales (R c 3> r c ), which allows to introduce the 
notion of quasicondensate characterized by the value no, 
the amplitude of the one-particle density matrix p(r) at 
intermediate distances r c -C r <C R c |§ , where 

p(r) = <*t(r,0)¥(O,0)>. (2) 

Local properties of the quasicondensate are identical to 
those of the genuine condensate, which justifies the 1 /mi- 
effect. However, in contrast to the 3D case, non-zero no 



in two dimensions is only due to finite interparticle in- 
teraction, which, in its turn, reduces the drop of K m . 
Another complication is associated with the fact that 
fluctuation region near T c (or n c if temperature is kept 
fixed) is rather wide. Also, the gas parameter is usually 
not so small in 2D. Thus, one may observe rather broad 
transformation of local correlation functions near T c or 
n c . 

The l/m!-effect itself may be used for the detection 
and study of BEC. So far experimental attempts to create 
2D waveguides for ultra-cold neutral atoms in specially 
arranged laser and/or magnetic fields (see, e.g., Refs. 
and Refs. therein) did not succeed in getting necessary 
conditions for BEC. Another promising system is spin- 
polarized hydrogen on helium film |6|,|^]. Small binding 
energy (~ 1 K), and strong derealization perpendicular 
to the film allow to consider the in-plane motion of hy- 
drogen as free. Recently it was announced || that the 
system undergoes BEC phase transition as the surface 
density is increased; the conclusion being based on the 
significant drop of the three-particle recombination rate. 
The advantage of such measurements in comparison with 
the standard search for the K-T transition in torsion-type 
experiments || is that they are not influenced by the sub- 
strate and are quasistatic. 

For realistic interatomic potentials and particle densi- 
ties the fluctuation region is too wide to allow adequate 
analytic treatment of correlation functions. We thus at- 
tempted quantum Monte Carlo (MC) simulation of the 
2D Bose gas in the grand canonical ensemble, by varying 
the degeneracy parameter nj (mT) (where m is the atom 
mass) through the chemical potential fi. Note that exper- 
iments with spin-polarized hydrogen on helium surface 
are done in the same setup, since surface density is con- 
trolled by \i of the bulk bufer gas. Surprisingly enough, 
under the same conditions the theory of Ref. jllj pre- 
dicts the decrease of K 3 by a factor of 400! The results 
presented below unambiguously prove that the theory of 
Ref. [|ll| is in error. 

Our MC simulation is based on the recently developed 
continuous-time Worm algorithm |l0| , which is extremely 
effective in calculating Green functions (at any tempera- 
ture), and is free of systematic errors. Thus apart from 
local correlation functions Eq. ([!]), we also evaluate the 
one-particle density matrix p(r) Eq. (||). 

We start with the definition of the model Hamiltonian 
on the square lattice 
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H = -t 



<ij> 



at a,4 



U 



= a\a,i 



(3) 



where a \ creates a boson on site i, and < ij > stands 
for the pairs of nearest neighbor sites. The particular 
form of the short-range interaction is not important in 
the dilute limit, and we restrict ourselves to the on-site 
Hubbard repulsion U. To ensure that underlying lattice 
plays no role, we choose parameters so that characteristic 
one-particle energies are much less than the bandwidth 
W = 8t, i.e., we require T,U <C 8t. Though the variable 
Yij = Yi—Tj is essentially discrete, in the quasi-continuous 
case, the density matrix p(r) is a smooth function of 
r =1 r l3 |. 

At the band bottom we may employ quadratic ex- 
pansion 6^ = k 2 /2m for the dispersion law with m = 
ft 2 /(2ta 2 ) — 1/2 (in what follows we use units such that 
% = 1, t = 1 and a = 1). The strength of interac- 
tion is naturally characterized by relative depletion of the 
condensate density at zero temperature £ = (n — no)/n, 
which can be easily derived in the framework of the stan- 
dard Bogoliubov transformation: 



£ = £//8tt 



(4) 



On one hand, we will study the case when £ is small. On 
another hand, the parameter U will be strong enough to 
see the effect of interparticle interaction on local correla- 
tors K m . 

To estimate the critical density we employ the uni- 
versal relation for the K-T temperature ijjjql , T c — 
irrig/2m — nnZ, where is the superfluid density at 
T ^ T r 



0. At a fixed temperature T the critical value 



of n s is given by 



= T/tt 



(5) 



For T = 0.2 and T = 0.1 we have nj « 0.064 and nj w 
0.032, correspondingly. Though the critical density n c > 
rig, the difference is normally not large. Hence, for an 
estimate of n c one can use rig (||). Note that for both 
above-mentioned temperatures (for which we will present 
below our numerical results) the K-T transition occurs at 
relatively low densities. 

A typical evolution of p(r) with increasing degener- 
acy is presented in Fig. 1. The curves are obtained 
for the lattice with L x L = 80 2 sites (we use peri- 
odic boundary conditions). In a normal state [Fig. 1(a)] 
there is only one characteristic length-scale - thermal de- 
Broglie wavelength, and p(r) decays exponentially with r. 
In Figs. l(c,d) we see a well-developed quasicondensate 
state, characterized by two different length-scales - after 
fast short-range decrease to a certain value (quasiconden- 
sate density no), p continues to decay very slowly. The 
case shown in Fig. 1(b) is an intermediate one. Though 
there is no pronounced bimodal shape yet, the decay of 



p at larger r's is anomalously slow. This type of behav- 
ior is observed in a rather wide range of variation of the 
degeneracy parameter (fluctuation region). 
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FIG. 1. Evolution of one-particle density matrix with den- 
sity. T — 0.2, U =1.0. Solid curves in Figs, (c) and (d) are 
obtained from Eqs. (^-^) , and contain no fitting parameters. 

Away from fluctuation region (at n > n c ) there exists 
an analytic representation for p(r) [fl3|| : Within Popov's 
functional-integral approach one obtains the follow- 
ing set of coupled equations that has to be solved self- 
consistently for p(r), no, and elementary excitation spec- 



trum E(k) = yje(k)[e(k) + 2n U], 

p(r) = e _A(r) j 5(r) , n Q = p(r -> oo) ; 
A(r) = 

p(r) -- 



[1 - cos(kr)] v{k) , 



(2tt 



rfk 



E{k) 
[1 — cos(kr)] 



(2n) 2 

e(fc) + n U - E{k) e{k)v{k) 



Here v{k) 



{ 2E(k) 
{exp[E(k)/T] - 



+ 



(6) 
(7) 

(8) 



U 



1 



E(k) J 
l} -1 is the Bose function, 

(9) 



U 



(mU/2n) ln(l/d fc c ) 



is the effective interaction, do is a cutoff for distance, and 
k c is a typical momentum. Expression (|^) implies that 
dok c -C 1. When the degeneracy parameter is on the or- 
der of unity or larger, to a good approximation one may 
set in Eq. (0) k c ~ y/n. In our model, do is just the in- 
tersite distance. The function p(r) describes short-range 
decay of the density matrix to the quasicondensate den- 
sity value no- The long-range decay of p is described by 
slowly growing exponent A(r). For n > n c , we observe a 
remarkably good agreement between p(r) calculated from 



2 



Eqs. (|6j-|9|), and our MC results; the agreement becomes 
progressively better away from the fluctuation region [see 
Figs. l(c,d)]. 
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FIG. 2. The two-particle density correlator as a function of 
density at U = 0.4. 
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FIG. 3. The three-particle density correlator as a function 
of density at U = 0.4. 

We now turn to the local density correlators ([!]). In 
Figs. 2 and 3 we present the data for Ki and K3 as, 
functions of density n for two temperatures T = 0.2 and 
T = 0.1, at U = 0.4. [The system size ranges from 
80 2 , for higher n's, up to 300 2 , for lower n's.] We see a 
pronounced strong decrease of K m when density varies 
from n < n c to n > n c . The most striking result is the 
very broad cross-over region. Both correlators start to 
decrease at densities well below n c . 

It is important to verify that our results are not ar- 
tifacts of finite-size effects since for finite L there exists 
a considerable fraction of genuine condensate even at fi- 
nite temperature. (Almost by definition, this fraction is 



given by p(r*,), where ~ L.) Hence, the change in lo- 
cal correlators could be due to global condensation (like 
in 3D case), rather than quasicondensation. We checked 
explicitly that the point U = TO, n = 0.04, T = 0.2, 
where there is already a pronounced decrease of Ki and 
K3, is not sensitive, within the error bars, to the system 
size for L = 60, 100, and 200. Also, we proved that at 
this point the value of /o(r*) is not a relevant quantity, 
being very small at our largest available L's (see Fig. 4). 
Most convincingly, in Fig. 4 we see that the crossover in 
Ki starts well before p(r*) becomes appreciable. 
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FIG. 4. K 2 (n) and p{r,,n) curves for U = 1.0, T = 0.2, 
demonstrating decoupling of short-range and long-range cor- 
relation properties (here r* = 40). 

The decrease of K m in the region n < n c is indicative of 
strong quasicondensate fluctuations in the normal state. 
In principle, this behavior is not unexpected, since no is 
not directly related to the superfluid density ns- If the 
concentration of vortices is small in the fluctuation region 
n < n c , the quasicondensate can survive, being related to 
the short-range correlation properties. The variation of 
no thus has a form of cross-over rather than a transition. 

In the limit of very small (but finite) interaction K m 's 
should change their values from ml to 1 throughout the 
transition. The data presented in Figs. 2 and 3 demon- 
strate two characteristic plateaus at n < n c and n 3> n c , 
the ratio between the two (at T = 0.1) being equal to 
w 4.6 for K3, and to w 1.85 for K2, and smaller than 
to!. The absolute values of correlators are also consider- 
ably smaller than in the case of negligibly small U. We 
thus conclude that even U = 0.4 is not small enough to 
yeald an idealized picture. Fig. 5 clearly demonstrates 
decreasing amplitudes of the 1/2-effect with increasing 
interaction. A similar picture was also observed for the 
correlator K3. At this point we note that K(U) depen- 
dence is not universal and can be sensitive to a particular 
form of the interaction potential. It is crucial however 
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that we observe weakening of the l/m!-effect with in- 
creasing U, which is in a sharp contradiction with Ref. 
p[ predicting an enormous enhancement of the effect. 
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FIG. 5. K2(n) curves for various coupling strength U, cal- 
culated at T = 0.2. 

The three-body dipole recombination rate W3 |T^ 1 
(proportional to K3) of spin-polarized atomic hydrogen, 
adsorbed on the surface of superfluid helium was mea- 
sured recently in Ref. Hydrogen dynamics perpen- 
dicular to the surface is quantized since there is only one 
localized state, and at relevant temperatures (~ 200 mK) 
the system may be considered as purely two-dimensional. 
In these experiments W3 was studied as a function of sur- 
face density It was found that recombination rate 
falls drastically with increasing n, and the decrease of W3 
starts well before the critical point n c is reached. This re- 
markable result is in qualitative agreement with our MC 
simulations. Unfortunately, direct quantitative compari- 
son is not possible because hydrogen atom derealization 
perpendicular to the surface is density dependent. In- 
deed, n can not exceed some maximum value n m . When 
n — * n m the absorption energy goes to zero and the hy- 
drogen atom wave-function in the direction perpendic- 
ular to the surface essentially changes its form due to 
collective effects jl6). This specific restructuring of the 
wavefunction leads to the additional drop of W3 (through 
the density-dependent factor a(n) in W3 = aK 3 ). As a 
result, if n c is close to n m , the observed rate W3 may 
drop by a factor much larger than 6 (up to 40) |ll|. As 
far as we can see, this is the only possible explanation for 
the measured ratio W 3 (n <C n c )/W 3 (n »ri c ) >6 |. 

Summarizing, we studied quasicondensation in a two- 
dimensional interacting Bose system in a rather interest- 
ing and experimentally important regime when interac- 
tion is not very small. We traced the evolution of one- 
particle density matrix and local correlators with increas- 
ing degeneracy parameter. We found that quasiconden- 
sate features appear far away from the K-T transition 



point. The effect of qusicondensation on local correla- 
tion properties - the main phenomenon directly relevant 
to the experiment - is clearly seen in this region, but its 
strength is rather sensitive to the interparticle interac- 
tion. 
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